Dipole and Bloch oscillations of cold atoms in a parabolic lattice 
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The paper studies the dynamics of a Bose-Einstein condensate loaded into a ID parabolic optical 
lattice, and excited by a sudden shift of the lattice center. Depending on the magnitude of the 
initial shift, the condensate undergoes either dipole or Bloch oscillations. The effects of dephasing 
and of atom-atom interactions on these oscillations are discussed. 



1. Bloch oscillations (BO) of a quantum particle in a 
periodic potential are one of the most fascinating phe- 
nomena of quantum physics Q. Since the pioneering 
experiment |2| in 1996, this phenomenon has been inten- 
sively studied for cold atoms in optical lattices @, with 
recent emphasis on quantum statistical (Fermi or Bosc) 
and atom-atom interaction effects. In particular, the dy- 
namics of degenerate Bose gases, on which we will focus 
here, was studied experimentally in [3- 0- E3l ■ ^ should be 
stressed from the very beginning that, when addressing 
this problem theoretically, one has to distinguish between 
quasi one-dimensional lattices (created by two counter- 
propagating laser beams) and truly ID lattices (or so- 
called modulated quantum tubes). Indeed, in the former 
case the number of atoms per well of the optical lattice 
can be as large as 10 3 — 10 4 , and a mean field approach 
(based on the Gross-Pitaevskii or nonlinear Schrodinger 
equation) is generally justified. This is not the case of the 
truly ID lattices, where only few atoms occupy a single 
well, and, hence, a microscopic analysis is required. For 
a tilted infinite lattice such analysis .based on the Bose- 
Hubbard model, was presented in 0, H, 0, where two 
regimes of BO - quasiperiodic and irreversible decaying 
- were identified. 

When referring to the typical laboratory experiments, 
an additional complication stems from the harmonic con- 
finement along the lattice. Clearly, harmonic confine- 
ment should modify BO of bosonic atoms, and the aim 
of this work is to estimate its effect. At the same time, 
parabolic lattices have their own interest, because they 
allow to study dipole oscillations of BECs. Recent ex- 
periments [10| have shown that there is a fundamental 
difference between dipole oscillations in quasi- and truly 
ID lattices. While in the former case the main effect of 
the periodic potential can be taken into account by sim- 
ply substituting the atomic mass by its effective mass in 
the ground Bloch band one observes a rapid decay 
of oscillations in the latter case. In the present paper we 
also briefly discuss dipole oscillations of a BEC in truly 
ID lattices, parti ally overlapping in this part with recent 
theoretical work |12| . 

2. We consider atoms in a parabolic lattice potential 
V(x) — Muj 2 x 2 /2 — V cos 2 (2irx/d), where the atoms are 
set into motion by a sudden shift of the trap origin. The 
relevant parameters of the system are the hopping matrix 
element J (defined by the amplitude of the periodic po- 
tential V), the 'parabolicity' v = Mio 2 d 2 , and the initial 
shift Iq — Ax/d. Using the single band approximation 



and neglecting atom-atom interactions, the dynamics of 
the system is described by the pendulum model [3J, 1 1 3| , 

iMi = -j^ai - ( a; +! + a '-i) > W 

where a/(i) is the complex amplitude of the atoms in 
the l-th well of the optical lattice. The separatrix of the 
pendulum corresponds to the shift 

r = 2{J/u) 1 / 2 . (2) 

If the initial shift Iq < 1*, the pendulum shows oscilla- 
tions around the equilibrium point and, referring to the 
original system, this regime is regarded as that of dipole 
oscillations of the atoms. If > I*, the pendulum is in 
the rotational regime, and the dynamics of the atoms can 
be regarded as BO in a local static field F = vl^jd. 

3. We begin with analyzing the regime of BO for an 
ideal bosonic system (i.e., no atom-atom interactions). 
Because the local static force F is not homogeneous, 
one has an additional process of dephasing of BO in a 
parabolic lattice, as compared to the paradigmatic case 
of a homogeneously tilted lattice. When discussing the 
mean atomic momentum, we can estimate this effect of 
dephasing by evaluating the sum 

p(t) ~ cxp(-m 2 /27 2 ) sin[(cj B + vm/Ti)t] , (3) 

m 

where u>b = vIq/Ti, m = I — Iq, and 7 is the width of 
the atomic wave packet (measured in units of the lattice 
period). Substituting the sum by an integral, we obtain 

p(t) ~ exp(-t 2 /2r 2 ) sm(uj B t) , (4) 

where 

r 7 = 71/71/ . (5) 

It is seen in Eq. that the dephasing time t 7 is de- 
fined by both the wave packet width and by the trap 
frequency (see Fig. 2] below). On the basis of this re- 
sult one might conclude that a narrow wave packet is 
preferable for studying BO in parabolic lattices. This is, 
however, not exactly true because a narrow wave packet 
implies a lower contrast of the interference pattern mea- 
sured in laboratory experiments. Thus, one has to keep 
a compromise between the contrast and dephasing, when 
preparing the initial wave packet. 
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FIG. 1: Bloch oscillations of N = 5 atoms in a parabolic 
lattice with parabolicity ^ = 0.04J: (a) mean momentum 
of non-interacting atoms; (b) mean momentum of interacting 
(W = 0.2 J) atoms; (c) macroscopic coherence of the interact- 
ing atoms. Initial shift lo — SI* = 80. 



The irreversible decay of BO according to Eq. Q is 
a consequence of our approximation of the sum by an 
integral. Without this approximation, the decay of os- 
cillations is followed by periodic revivals with a period 
T„ = 2nT%/v [Tj]- One of these revivals is illustrated in 
the upper panel of Fig. ^ which shows the dynamics of 
the mean momentum of the non-interacting atoms in the 
parabolic lattice with parabolicity v — 0.04J. As initial 
state of the system we choose here the ground state of 
the atoms in a parabolic lattice with a slightly tighter 
confinement v' — 4^, which was then shifted by a dis- 
tance lo = 81* = 80. Note that by changing v' we change 
only the dephasing time [through the change of the wave 
packet width 7 = 7(V)] , while the revival time is defined 
exclusively by the parameter v. 

4. Next we address the effect of atom-atom interac- 
tions. In the case of quasi one-dimensional parabolic 
lattices, Bloch and dipole oscillations of the interacting 
atoms were studied in a number of papers, using the 
mean- field approach El El El El El OS As known, 
the mean-field approach is justified in the limit of large 
occupation number n — > 00 and vanishing microscopic 
interaction constant W — > 0, and leads (in the simplest 
case) to the discrete nonlinear Schrodinger equation, 

thai = ^ 2 a z - — (a t+ i + a/_i) + g\m\ 2 ai , (6) 

where g — WN is the macroscopic interaction constant. 
In our present work, we focus on the case of truly one- 
dimensional lattices, where the mean occupation number 
n *~ 1. Clearly, the mean field approach is not applicable 
here and one has to treat the system microscopically by 



using, for example, the Bose-Hubbard model, 

1 \ 1 I 1 

The main question we address below is the effect of atom- 
atom interactions on the Bloch dynamics depicted in the 
upper panel of Fig. ^ 

First we shall discuss the initial conditions in some 
more detail. Throughout the paper we shall consider 
the ground many-body state of the atoms in a parabolic 
lattice as the initial wave packet (which is shifted then 
by the distance lo). Clearly, along with the ratio J/v 
this state is also defined by the ratio of the interaction 
constant to the hopping matrix element. Namely, it is es- 
sentially given by the symmetrized product of the single- 
particle atomic state for W < J , while it may resem- 
ble the Mott-insulator state for W 3> J 12]. In what 
follows we restrict ourselves by a relatively weak interac- 
tion. Then the ground state of the system can be well 
approximated by the many-body wave function 

where |n) = | . . . , rt_i, no, n>x, . . .) is the Fock basis and 
the ai satisfy the stationary nonlinear Schrodinger equa- 
tion 

7^l 2 ai ~ fa^ 1 + a l- 1 ) + 9\ a i\ 2a i = E oai ■ (9) 

For example, for N = 5, v = 0.04J and W = 0.2J, 
the overlap of the state (JSJ) with the exact ground state 
|*o) is |(*o |*o) I 2 = 0.97. We note that the state © is 
completely coherent and is analogous to the super-fluid 
state in a homogeneous lattice. We shall characterize the 
macroscopic coherence of the given many-body state | '&) 
by the maximal eigenvalue A of the single-particle density 
matrix 

pi, m = N- 1 (*\a+a m \y) . (10) 

Then the macroscopic coherence of the state (JSJ is A = 1 . 

We proceed with the dynamics. The middle panel in 
Fig. ^ shows the mean momentum of N = 5 interacting 
atoms (W = 0.2 J). In comparison with the noninteract- 
ing case (upper panel), a qualitative change is noticed. 
This change can be understood by analyzing the macro- 
scopic coherence of the system, shown in the lower panel. 
It is seen that the macroscopic coherence oscillates with 
some characteristic period T\y- In the case of a homo- 
geneously tilted lattice these oscillations were studied in 
Ref. 0- The origin of the oscillations was shown to be 
the Stark localization of the single-particle wave func- 
tions which, together with the discreetness of the atom 
number, leads to the following expression for the macro- 
scopic coherence, 

A = exp(-2n[l - cos(Wt/h)]) . (11) 
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FIG. 2: Dipole oscillations of N — 4 atoms in the parabolic 
lattice with parabolicity v = 0.04J: (a) mean momentum 
of non-interacting atoms; (b) mean momentum of interacting 
(W = 0.2J) atoms; (c) macroscopic coherence of interacting 
atoms. Initial shift lo = Z*/2 = 5. 

In Eq. (|llfl n is the mean number of atoms per lattice site 
pcj and the limit Fd ^ J is implicitly assumed. Since 
for the considered local static force Fd = vIq = 3.2J 
Stark localization is not complete, the oscillations of the 
macroscopic coherence decay in time. Nevertheless, if 
this irreversible decay of coherence is slow on the time 
scale of the dephasing time, one can observe the revival 
of BO of the interacting atoms - an effect which attracts 
much attention because it provides an independent and 
accurate method for measuring the microscopic interac- 
tion constant W. 

5. Let us now turn to the case Iq < I*. Here we 
meet dipole oscillations of a BEC with a characteristic 
frequency given by the frequency of small pendulum os- 
cillations ujq — {vJ) 1 ' 2 /%. (We recall in passing that the 
frequency of BO was given by lob = vIq/% ~ 2u>qIq/1*, 
Iq 3> I*.) For vanishing atom-atom interactions these 
dipole oscillations are shown in the upper panel of Fig. [2 
where lo = l*/2 = 5, and the other parameters are the 
same as in Fig. ^ The dephasing time t 7 is again given 
by Eq. JSJ but with the parameter v substituted by the 
nonlinearity parameter v = z//8 |2l|. (The latter param- 
eter also defines the revival time.) The middle and lower 
panels in Fig. [21 refer to interacting atoms. An exponen- 
tial decay of the macroscopic coherence is noticed. The 
other point to which we want to draw the attention of 
the reader is that a moderate interaction stabilizes the 
dipole oscillations against dephasing. Within the mean- 



field approach (which reduces the Bose-Hubbard model 
to the discrete nonlinear Schrodinger equation) , this phe- 
nomenon is discussed in Ref. [lflj . 

6. In conclusion, we have shown that the dynamics 
of cold atoms in parabolic lattices is governed by the 
relation between two characteristic times - the dephasing 
time r 7 and the decoherence time Tw 

The dephasing time is inversely proportional to the 
width 7 of the initial wave packet and the nonlinearity v, 
which, in turn, is defined by the initial shift Iq of the wave 
packet relative to the separatrix I* = 2(J/v) 1 / 2 . Namely, 
v = v/8 for Iq <C I*, and v = v for Iq 3> I*. It is interest- 
ing to estimate the dephasing time in a typical laboratory 
experiment. Taking, as an example, the recent exper- 
iment [lfj with rubidium atoms in an array of axially 
modulated quantum tubes, we have v = 0.0014-Er and 
J = 0.38Er for the modulation amplitude (depth of the 
optical lattice) of one recoil energy. This gives a separa- 
trix position I* = 33, and a period To = 12.1 ms of small 
dipole oscillations. Assuming a dilute gas (which, in fact, 
is not the case realized in the cited experiment) the width 
of the initial wave packet is 7 ~ (J/4^) 1 / 4 = ^/J* '/2 » 3 
and, hence, the dephasing time r 7 = 85 ms for dipole 
oscillations, and r 7 = 10.6 ms for BO. Note that these 
are upper estimates for the dephasing times, and for the 
initial shift Iq closer to the separtrix the dephasing times 
are essentially smaller. It is also worth noting that there 
is a maximal shift Iq above which the single band ap- 
proximation (used throughout the paper) is not valid. 
The crucial parameter here is the energy gap between 
the Bloch bands (A = Q.§Er for the specified parame- 
ters). The analysis of BO in a parabolic lattice beyond 
the single-band approximation will be subject of a sepa- 
rate paper. 

The decoherence time tw is defined by the character- 
istic density of the atomic gas n and by the value of the 
microscopic interaction constant W. The latter, in turn, 
is defined by the s-wave scattering length and by the 
degree of confinement of the atoms in the wells of the op- 
tical potential. In particular, in the experiment |T3, the 
quantum tubes were created by two crossing quasi ID 
optical lattices with an amplitude V = 30-Er. For the 
axial modulation with V = this gives W — 0.73Er. 
For this relatively high value of the interaction con- 
stant few atoms per one tube are enough to destroy the 
dipole/Bloch oscillations on a very short time scale. This 
qualitatively explains the results of the experiment |10| . 
where the number of atoms per one quantum tube was 
around 20. To observe the effects discussed in this paper 
one has to decrease either the atomic density or the in- 
teraction constant (e.g. by use of a Fishbach resonance), 
as compared to those of Ref. • 
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